Majorana fermions in carbon nanotubes 
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We show that carbon nanotubes (CNT) are good candidates for reaUzing one-dimensional topo- 
logical superconductivity with Majorana fermions localized near the end points. The physics behind 
topological superconductivity in CNT is novel and is mediated by a recently reported curvature- 
induced spin-orbit coupling which itself has a topological origin. In addition to the spin-orbit 
coupling, an important new requirement for a robust topological state is broken chirality symme- 
try about the nanotube axis. We use topological arguments to show that, for recently realized 
strengths of spin-orbit coupling and broken chirality symmetry, a robust topological gap ~ 500 mK 
is achievable in carbon nanotubes. 
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Introduction: Topological superconductors (TS) in 
dimensions £' = 2,1 with broken time-reversal (TR) sym- 
metry have recently attracted a lot of attention [iHlOj. 
These systems support Majorana fermion excitations at 
order parameter defects such as vortices and sample 
edges. Majorana fermion excitations, with second quan- 
tized operators 7 satisfying the self-hermitian condition 
7^ = 7, can be construed as quantum particles which are 
their own antiparticles fTl| . The self-hermiticity of Majo- 
rana fermions (MF) leads to a 2D quasiparticle exchange 
statistics which is non- Abelian [l^] . The non- Abelian 
statistics of MP's can be used as a robust quantum me- 
chanical resource to implement fault-tolerant topological 
quantum computation (TQC) l3, 14 1. 

In ID TS with broken TR symmetry, Majorana 
fermion modes are supposed to be trapped at the two 
ends of a quantum wire which, in a 2D quan- 

tum wire network [l5| . can potentially lead to success- 
ful demonstration of non-Abelian statistics as well as 
TQC [H-li^. Recently, a number of systems with ef- 
fective Rashba-type spin-orbit coupling and broken TR- 
invariance have been proposed to be in an appropriate 
TS state with localized MF modes In this paper 

we add to this list an attractive new candidate - carbon 
nanotubes (CNT) - which have already become impor- 
tant experimental systems due to their remarkable elec- 
tronic and transport properties [l8l - [20| . Note that the 
superconducting proximity effect from an s-wave super- 
conductor on carbon nanotubes has also been realized 
[? ]. By employing robust topological arguments we 
show that CNTs in proximity to ordinary superconduc- 
tors such as Al, Nb present a good platform for realizing 
Majorana fermions at the ends of the nanotube provided 
the chirality symmetry about the tube axis can be bro- 
ken. Using parameters from recent experiments 3, 2^ 
we show that, under appropriate external conditions, it 
is possible to realize a TS state gap protecting the Ma- 
jorana fermion excitations as high as ^ 500 mK. 

The physics behind the topological state in CNT is 
novel and is mediated by a spin-orbit coupling (SOC) 



which itself has a topological origin. SOC in carbon- 
based systems (such as Graphene) is generally negli- 
gible due to a small ('^ 8 meV) atomic SOC of car- 
bon atoms which gets further suppressed by the band- 
structure. Nonetheless, CNT, which is cylindrically 
roUed-up Graphene, can have a reasonably strong SOC 
due to curvature- induced topological effects [19r.22.] . The 
band structure of CNT inherits the two momentum space 
Dirac cones of Graphene, but the momenta along the cir- 
cumferential direction are discrete and quantized due to 
finite size effects. Electronic states with quantized cir- 
cumferential momenta arising from different Dirac cones 
have opposite chirality (indexed by chirality index = 
1,-1 in Eq. ([1])), which indicates the clockwise or an- 
ticlockwise circulation about the nanotube axis. In one 
circulation around the circumference the electronic or- 
bitals undergo a complete rotation about spin, which re- 
sults in an effective topological spin-orbit coupling with 
strength inversely proportional to the tube diameter d 

MB- 

It has been shown recently that for tubes with 
diameter d = 5 nm, the spin-orbit energy splitting Ago 
can be as high as ~ 0.4 meV. Note that A50 arises from a 
topological Berry phase accumulated by orbital rotation 
around the circumference and has no explicit momentum 
dependence as is customary in Rashba or Dresselhaus 
couplings in ordinary semiconductors. We show below 
that this type of spin-orbit coupling, even in the absence 
of an explicit momentum dependence, can produce a one- 
dimensional TS state provided time-reversal and chiral 
symmetries are also broken in CNT. 

The topological spin-orbit coupling and broken TR in- 
variance (due to an external transverse magnetic field) 
are not sufficient by themselves to realize a finite gap in 
the TS state in the nanotube. One more crucial ingredi- 
ent for a non-zero TS state gap is broken chiral symmetry 
about the tube axis. In the presence of chiral symmetry 
(i.e., rotational symmetry about the tube axis) all states 
in the nanotube can be classified as circulating clockwise 
or anticlockwise about the tube axis (i.e., CTz in Eq. ([1]) 
is a good quantum number). In one of the main results 
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of this paper we show on quite general grounds that in 
this case the gap in the topological state, even if the 
topological invariant is non-trivial, identically vanishes. 
Therefore, for a finite TS state gap with robust Majo- 
rana fermions, the nanotube Hamiltonian must mix and 
split the chiral eigenstates. Such couplings between the 
clockwise and the anticlockwise states can be induced by 
breaking the rotational symmetry about the nanotube 
axis using the superconducting proximity effect itself or it 
can arise from disorder (which can also break this symme- 
try). In a recent paper Jespersen et al. has shown ^22] the 
existence of such a term /S.k.k' in a multi-electron CNT 
with a magnitude of ^k,K' of the same order as A50. 
By calculating the ap pro priate topological invariant (the 
Pfaffian invariant 23j,|24|) and the bulk quasiparticle gap 
we show that, in such a nanotube, a finite TS state gap 
~ 500 mK is achievable with a modest value 1 T of 
an external magnetic field. Our calculations and results 
establish the single walled carbon nanotube as an attrac- 
tive new candidate for topological physics and Majorana 
fermions potentially leading to TQC. 

Hamiltonian for carbon nanotube and topological spin- 
orbit coupling: The states in a carbon nanotube are 
characterized by a chirality index cr^ ±1, a pseudo-spin 
Pz and a parallel spin component Sz. The Hamiltonian 
of a CNT in the presence of a transverse magnetic field 
B±, but no chirality breaking term, is given by, 

HcNT = VFcrz{{k\\+Kiiaz)p.j; + {ko + aazSz)py} + Bj_s^ 

(1) 

where the Hamiltonian, for B± = 0, commutes with the 
anti-unitary time-reversal operator T = iSyUxK which 
satisfies = —1. Here Px,y are Pauli matrices in the 
sub-lattice degree of freedom, while cr^ is the Pauli ma- 
trix whose eigenstates arise from different valleys isT, K' 
of the underlying Graphene sheet. The wave- vector /c| 
is the wave-vector of the electron along the axis of the 
carbon nanotube. Since the Dirac cones if, K' of the un- 
derlying Graphene sheet occur at wave-vectors with op- 
posite signs, they have opposite projection K^^Gz — 
along the axis of the nanotube. Confinement along the 
circumference of the nanotube leads to a minimum cir- 
cumferential momentum of magnitude kg and a corre- 
sponding minimum gap vpko in the Dirac cone where 
vp is the fermi velocity of the underlying Graphene 
sheet. Electrons near the bottom of each band near K 
or K' have circumferential group velocities in the clock- 
wise or in the anti-clockwise direction (corresponding to 
az = ±1) respectively. The atomic spin-orbit coupling of 
carbon leads to electrons with different spin orientations 
to experience different Aharonov-Bohm fluxes, shifting 
ko ^ ko + aSzO'z, where a is the topological spin-orbit 
coupling strength [l^. Finally B± is the Zeeman split- 
ting of the electrons from an applied magnetic field in the 
transverse (x-) direction. For the moment we are ignor- 
ing any chirality breaking term in the Hamiltonian (so cr^ 



is a good quantum number). This is because we want to 
prove below one of the central results of this paper which 
states that, in the presence of chirality symmetry, there 
cannot be a non-zero TS state gap even when the Pfaf- 
fian topological invariant given by Eq. is non-trivial. 
It follows that, for a robust TS state in CNT, mixing of 
the clockwise and the anti-clockwise states in the nan- 
otube is a crucial requirement as is a non-zero spin-orbit 
coupling a. 

BdG Hamiltonian and the Pfaffian invariant: In the 
presence of an s-wave superconductor proximity coupled 
to the nanotube, a superconducting pair potential A will 
be induced on the CNT. In this case, the singlet BdG 
Hamiltonian for the CNT that anti-commutes with the 
particle-hole operator A — —iTyT is given by. 



H 



BdG 



VF{k\\crzPx + K\\px + koazPy + aSzPy) ~ p 
B^Sx + Atx. 



Tz 

(2) 



g = sgn P/(Ai7sdG(fc|| 



The BdG Hamiltonian for the CNT is topologically non- 
trivial and supports Majorana fermions at the ends pro- 
vided the Pfaffian invariant written as, 

0)))x(p/(Ai/BdG(A:|| 

(3) 

is — 1 2j]. However, in addition to the non-trivial value of 
Q, for a useful topological state with localized Majorana 
fermions protected from thermal decoherence, the CNT 
must have a non-zero (proximity-induced) topological su- 
perconducting gap. Below we show that, for Q = —I, a 
non-zero TS state gap is not possible as long as the chi- 
rality (given by the tr^ in Eq. ([l])) is a good quantum 
number. 

Absence of topological gap in the presence of chiral 
symmetry: The BdG Hamiltonian in Eq. ([2]) commutes 
with az so that different chirality sectors az — H are 
decoupled. Here we show that such Hamiltonians are 
necessarily gapless whenever it has a non-trivial Pfaffian 
invariant. To do this, we first assume that Q = — 1 so 
that the Pfaffians of Hsdo at the particle-hole symmet- 
ric points fc|| = 0;7r have opposite signs. Next construct 
another BdG Hamiltonian, 

Hsdc = vpikiipx + K\\px + koazPy + aSzPy) ~ p Tz 
+ B^Sx + Atx, (4) 

such that Hsdcikw) is particle-hole symmetric at each 
value k\\ (i.e., not just at fc|| = 0, ±tt as in Eq. ([2])). This 
allows us to define a Pfaffian for HsdC at each value of 
k\\ and not just at /c|| — 0, ±7r. Note also that H[BdG{k\\) 
has the same spectrum as Hsdcik^^) and, in addition, 

-H^BdG(fe||=o) = HBdG(k^^=o) and the equahty also holds 
for k\\ — ±7r. Therefore, if Q = — 1 for HBdG, so is Q for 
Hsdo- Since we have assumed that HBdG does indeed 
have Q = —1, it follows that the Pfaffians of HsdG at 
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FIG. 1: Bandstmcture of a (9,3) chiral CNT. Such CNTs have 
a band minimum away from fc|| = and therefore the bands 
have a linear intersection near fe|| =0. 

fc|| — 0,TT also have opposite signs. Therefore the Pfaffian 
of HbcIG must vanish at some fc|| between fc|| = and 
fc|| = TT. Note that since Hsdc is particle-hole symmetric 
at all points fc|| between and tt (but HbcLG is not) we 
can define a Pfaffian for HbcIG at arbitrary fc|| and it 
must be zero somewhere on the fc|| line for Q = —1. The 
Pfaffian of a BdG Hamiltonian can vanish at a fc| | point 
only when one of the eigenvalues vanish at that point. 
Further, since commutes with both HsdG and HBdGi 
all eigenstates have definite values of ctz- Therefore, if 
HBdG is gapless at some value of k\\ and Gz, so is HBdG- 
This shows that carbon nanotubes are gapless at some 
value of fc|| whenever the Pfaffian invariant is non-trivial, 
Q = -1. 

Chiral symmetry breaking and nanotuhe hand- 
structure: To obtain a non-zero topological gap when 
the Pfaffian invariant Q = —1 one must add a chirality 
breaking term to the nanotube Hamiltonian which does 
not commute with cr^. Such a term can arise from 
disorder effects and/or broken rotational symmetry 
about the nanotube axis due to the substrate. Note that 
the superconducting proximity effect itself can be used 
to break the rotational symmetry about the CNT axis. 
With broken rotational symmetry about the tube axis 
the clockwise and the anticlockwise states are no longer 
eigenstates of the one-electron Hamiltonian and their 
splitting can be represented by the term IS.k,k' ■ Such 
a term, with a magnitude comparable to the energy 
splitting due to the topological spin-orbit coupling, has 
recently been experimentally observed Adding this 
term to Eq. ([1]) the nanotube Hamiltonian changes to, 

HcNT = VF(Jz{{k\\+ K\\az)px + {ko + a(JzSz)py} + B^s.^ 
+ l^K,K'<yx- (5) 

To understand the effects of ^k,k' let us first ignore 
its contribution to the band structure. The dispersion of 
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FIG. 2: The combination of t^KK' opens up a gap in the 
dashed band-structure with an odd number of fermi surfaces. 
We have chosen a = 400/ieV,-B^ = 400/ieV (corresponding to 
0.5 T) and A^k' = 200^ieV. 

the conducting states in the nanotube are then given by 



e(fc, cr, s) = VF\J{k\\ - K\\(JzY + (fco + ascjzY - p, (6) 

where we have used the fact that the Hamiltonian com- 
mutes with az- For simplicity, we consider a nanotube 
in the semiconducting regime so that we can focus only 
on the states in the conduction band. The correspond- 
ing band-structure near the relevant intersection point at 
fc = is plotted in Fig. [T] Each curve in Fig. [1] actually 
represents two bands (because of near spin-degeneracy) 
split by much smaller energy scales B± and a, and these 
are shown as dashed lines in Fig. [2] 

From the dashed lines in Fig. [2] it is clear that even 
though the conduction bands {E > 0) are split by the 
topological spin-orbit coupling a and the transverse field 
B±, all values of the chemical potential give rise to an 
even number of Fermi surfaces. We know, on the other 
hand, that for a topological superconductor with an odd 
number of Majorana fermions at each end of the nan- 
otube it is necessary to obtain an odd number of bands 
at /c|| =0 [1^. The solid curves in Fig. [2] depict the 
band-structure near fc|| = with a non-zero chirality 
breaking term Ak.k' in the nanotube Hamiltonian. The 
term Ak,k' opens up a gap in the one-electron spectrum, 
but, in contrast to the Rashba-coupled systems, chemi- 
cal potential in the gap at =0 does not produce an 
odd number of Fermi surfaces. There are in fact two 
ranges in (for the parameters in Fig. [2] on both sides of 
E 102 meV), for which the nanotube has an odd num- 
ber of Fermi surfaces and, consequently, is in a gapped 
topological state in the presence of a proximity-induced 
pair potential A. 

Topological state and Majorana fermions with chiral 
symmetry breaking: The BdG Hamiltonian of the nan- 
otube including the topological spin-orbit coupling a, 
transverse field B±, and the chirality breaking term 
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CNTs with end-state Majorana fermions, the chiral sym- 
metry about the tube axis must be broken. This can 
be achieved by the superconducting proximity effect it- 
self or by disorder which may also break the rotational 
symmetry about the nanotube axis. Using recently re- 
ported values of Ago and the chiral symmetry breaking 
term IS.k^k' we show that in carbon nanotubes a robust 
topological gap ^ 500 mK is achievable. 

J.S. thanks the Harvard Quantum Optics Center for 
support. S.T. thanks DARPA and NSF for support. 
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FIG. 3: Topological superconducting gap in carbon nanotube 
as a function of chemical potential /i. The gap is multiplied 
by the negative of the topological invariant Q (Eq. [3]), so it 
is positive in the nanotube TS state with Majorana fermions 
{Q = —1) and negative in the topologically trivial {Q = 1) 
superconducting state. We have chosen a proximity induced 
s-wave pair potential A — 300/ieV. 



^KK' can be written as, 



BdG 



VF{k\\crzPx + K\\px + koa^py + as^Py) 



^K,K'<yx — 



B I 



(7) 



For the Hamiltonian in Eq. [7] we first calculate the Z2 
invariant Q given in Eq. [31 In Fig. |31 we plot the su- 
perconducting gap (calculated by diagonalizing the BdG 
Hamiltonian in Eq. [T]) multiplied by the negative of the 
Z2 invariant Q as a function of the chemical potential 
p. Since Q = —1(1) signals the topologically non-trivial 
(trivial) superconducting state, positive values of this 
product indicate the TS state in Fig. |3l It is clear that 
there are well defined ranges of the chemical potential in 
which the nanotube supports Majorana fermions with a 
non-zero TS state gap ^ 500 mK. 

Conclusion: In this paper we show that carbon nan- 
otubes are an attractive candidate for realizing one- 
dimensional time-reversal breaking topological supercon- 
ducting state with zero-energy Majorana fermions local- 
ized near the ends. The physics of the TS state in CNTs 
is mediated by a novel, curvature-induced, spin-orbit cou- 
pling which is topological in itself. We show that, despite 
the presence of this coupling (reported in recent experi- 
ments 13, 2^ to be substantial, A50 ^ 0.4 meV), there 
cannot be a non-zero TS state gap in the nanotube in 
the presence of chirality symmetry even in a parameter 
regime where the relevant topological invariant is non- 
trivial. Thus, to obtain a robust topological state in 
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